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ABSTRACT 



The purpose of this project is to localize buried 

objects in water-saturated sand by variable incidence 

acoustic pulse reflections. In particular, thin aluminum 

plates are used to model the buried objects. The sound 

speed and the depth of the overlying sand layer are predicted 
2 2 

by using the T -X and the ray parameter methods. The 

experimentally computed results have an average 3.5% depth 

error. The calculated speed of sound in the water-saturated 
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sand layer differed by 8.1% between the T -X and the ray 
parameter methods . 
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INTRODUCTION 



I . 



A. OBJECTIVES 

Identification of buried objects in the ocean bottom is 
an important problem in many fields such as oil exploration, 
geophysical research, and ocean engineering. A simpler 
problem is the identification of layers of infinite extent 
and unknown properties in the ocean sediments. The objec- 
tives of this thesis are: 

1. to determine experimentally the depth of burial 

of a thin plate in a simulated sedimentary bottom; 

2. to determine the sound speed and the thickness 
of the overlying sand layer; and, 

3. to determine the sound speed and thickness of 
the thin plate. 

It was the intent of this thesis to determine all of these 
properties by remote sensing with acoustic pulses in order 
to simulate an ocean measurement with a ship and near- 
surface sensors. 

The first objective could have been accomplished by 
normal-incidence time-of-f light measurements; however, the 
second objective requires more information than that so a 
variable angle of incidence technique was used. This second 
objective was actually the principal goal since the third 
objective could not be achieved by using the technique 
described herein. The technique failed because the pulse 
could not be made short enough to distinguish the reflection 
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from the plate's lower surface. This is primarily the 
result of using very thin plates with respect to the acoustic 
wavelength. Perhaps the amplitude of the reflection from 
the top of the plate could be used to assist in identifying 
the plate material but the overlying sand layer makes this 
very difficult. 



B . APPROACH 

The approach used is as follows: 

1. Short pulses consisting of several cycles of a 
150kHz sinusoidal wave were used as the source 
pulse . 

2. A directional projector was used to minimize inter- 
ference from surface and side-wall reflections. 

3. The receiving hydrophone was mounted on a carriage 
so that the horizontal distance from source to 
receiver could be easily varied. 

4. Each of the thin plates were buried in the sand 
bottom of a small water tank. 

5. The plates were aluminum ranging in thickness from 
0.1 to 1.8 cm and were buried at depths from 1.8 
to 5.0 cm. 

6. Two analytical methods were used for reducing the 
data. One method was based on the relationship 
between the square of the travel time (T) of the 
reflections and the square of the source-to-receiver 
horizontal range (X) . The other method was based 

on maintaining a constant ray parameter (P) defined 
as the sine of the incident angle divided by the 
local sound speed. 



C. BACKGROUND 

Reference 1 gives an excellent historical development 
of various seismological methods for determining properties 
of sediment layers. The following is a brief summary of 
that account. 
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Karcher and Haseman made the first successful seismic 



reflection measurements in 1921. Techniques for acoustic 

determination of sediment properties have advanced greatly 

since then and generally fall into one of three categories: 

normal incidence, wide-angle reflection, and refraction 

methods. Most of the development has been done by the 

petroleum exploration industry; however, the geophysical 

community has used these methods to study the earth's crust 

as early as the late 1930 's. Dix, Krey and Durbaum 

refined the sound speed measurements for sediment layers 

below the first layer by properly accounting for the shape 

of the arrival time as a function of range curves . By 

1955, Dix had demonstrated the equivalence of T-X methc 
2 2 

and T -X method and, in 1968, Le Pichon refined the 
2 2 

T -X method until it became a standard analysis tool. 

A major enhancement was introduced by Bryan (see Ref. 2) in 

2 2 

1974 with the ray parameter method. Both the T -X method 
and the ray parameter method were used in the data analysis 
for this thesis. 
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II. THEORY 



A. T-X AND T 2 -X 2 METHODS 

A simple model of a layered sediment is shown in Fig. 
2.1. The layers have thicknesses H^ , and H 3 and sound 
speeds and . The T 2 -X 2 method [Ref. 1] can be 

applied to this situation as described in this section. 




Figure 2.1. Geometry of Layers 

For convenience, let the sound speed in the uppermost 
layer (the water layer) be less than that in the first 
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sediment layer. The length of the path that reflects from 
the top of the sediment layer is given by, 




X 2 + 4 (H 1 - s ) 2 



( 1 ) 



To convert this to travel time squared, divide through 
by V 2 : 



T 2 (X) = T 2 (0) + (— 4 ) X 2 

V 1 



( 2 ) 



where T(X) is the reflection travel time at horizontal range 

X and T(0) is the intercept of the travel time curve at 

2 

X = 0. Since T (0) is a constant, the reflection curve is 

a hyperbola. The corresponding curve for the direct path 

from the source to the receiver is a straight line and this 

hyperbola is asymptotic to that direct-path line. 

2 2 

When plotted in the T -X plane, however, these lines 

are both straight. In this way, the slope of the reflection 

curves can be easily determined. As can be seen from the 

2 2 

equation for the reflection path, the slope in the T -X 

plane is -4 which, in this thesis work, is calculated by 
V 1 2 2 

linear regression of the T -X data. This rather simple 
but accurate method is due to Green [Ref. 1] and is asso- 
ciated with the wide-angle seismic reflection technique. 

The paths that pass through one or more sediment layers 
before reflection (see, for example, R-^ i n Fig. 2.1) do not 
lead exactly to hyperbolas but the departure is not great 
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Figure 2.2. 



T-X and T 2 -X 2 



Curves 
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and only a small amount of error is introduced by fitting 

2 2 

this data with straight lines in the T -X plane. 

If the layers were extremely thick or had great differ- 
ences in sound speed, then the error might become apprecia- 
ble. The reciprocal of the slope of the R^ 2 curve is equal 
to the average speed (squared) through both the water and 
the first sediment layer: 

V = (H 1 +H 2 )/( (Hj/Vj.) + (h 2 /v 2 )) 

= (V 1 T 1 +V 2 T 2 )/(T 1 +T 2 ) (3) 

Representative T-X curves for this case are shown in 

Fig. 2. 2. a. Curve D represents the direct path and R^ 

represents the reflection curve for reflection from the top 

of the sediment. For R^ 2 , the exact curve (solid line) is 

given as well as the hyperbola approximation (dotted curve) 

2 2 

to that curve. The equivalent curves in T -X space are 

shown in Fig. 2.2.b. Here the inverse slopes give the sound 

speed squared. The intercept T(0) gives the two-way 

vertical travel time through the layers. 

2 2 

The principal limitation of the T -X method is shown 
in these figures. For any reflection other than the reflec- 
tion from the top of the sediment, the curve is not 

2 2 

exactly a straight line in T -X coordinates but it falls 
slightly below the equivalent straight line as range 
increases . 
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B. THE RAY PARAMETER METHOD 



In order to measure the properties of a layer without 
knowing the structure of the overlying layers, the ray 
parameter method can be used [Ref. 2J . Any ray path is 
specified by a fixed value of the ray parameter P defined 
by. 



P = sin 6 (z)/c (z) 



(4) 



where z is depth, c(z) is the sound speed profile and 0(z) 
is the angle that the ray makes with the vertical at depth 
z. Figure 2.3 shows the relationship between the incident 
angle and the ray parameter. 

According to the ray parameter method, the increment of 
travel time along the ray path gives the corresponding 
increments in depth and range: 



dz = c dTcos 0 



2 2 

(1 -Pc 



1/2 

c dT 



(5) 



dX = c dTcos 0 



P c 2 dt 



2 2 -V 2 

Pc (1 - P c ) dz 



( 6 ) 



After integrating over depth from the surface to depth 
z i , these expressions can be rewritten as follows. 



T (P) 



Z 1 

t - 1 ,, 2 2 , - 1/2 , 

J c (1 - P c ) dz 

0 



(7) 
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X 




X(P) 




2 2 

Pc(l -P c 



- 1/2 

dz 



( 8 ) 



These equations give the travel time T and the hori- 
zontal distance X along the ray path given by P in terms of 
integrals over depth. The equations must be solved for the 
sound speed variation c(z). This can be done for some 
special cases of layer structure. 

Consider a homogeneous layer (sound speed c^) extending 
from z^ to z^ +Az (Fig. 2.3). For the ray of a particular 
ray parameter value, there are two (in this case) rays that 
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return from the bottom with that same parameter value 
although they arrive at the receiver at different ranges. 
Note that these rays are parallel through the water column. 
The travel time and horizontal distance for these two rays 
are, according to Eqs . (7) and (8), 



AT (P) 



z +Az 

' - 1 ,, 2 2 , 

c (1 - P c ) dz 



AX ( P) 



2 



z„ +Az 

/ Pc (1 - P 2 c 2 



- 1/2 

dz 



(9) 



( 10 ) 



where the factor 2 accounts for the round trip of the ray. 
If the sound speed is constant at c^ in the layer, these 
equations can be integrated to give, 



AT 



AX 



-1 2 2 
2c 1 J '(l -Pc) 

2Pc^ ( 1 - P 2 c 2 ) 



- 1/2 

Az 



- 1/2 

Az 



(ID 

( 12 ) 



The expressions for sound speed and layer thickness are 
then , 



,1 AX, 
C 1 ( P AT 



1/2 



Az = i c 1 AT ( 1 -P 2 c 2 ) 



1/2 



(13) 



(14) 
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Figure 2.4. T-X Curve and the Ray Parameter 

The important feature of the ray parameter method is that 
the structure above the layer has dropped out because of 
the choice of ray path. 

In order to apply this technique in practice, an initial 
value of ray parameter was chosen based on a geometry that 
gives good clean reflection signals. After the T-X curves 
are drawn, this value of ray parameter, which gives the 
slope of the T-X curves at the points of interest, was used 
to identify the two parallel return rays. The principal 
difficulty of this method is that the T-X curves are not 
straight lines and, therefore, curves must be fit to the 
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experimental data. This introduces some uncertainty into 
the process but, fortunately, the values for layer thickness 
is not very sensitive to the type of fit that is done. 
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III. EXPERIMENTAL PROCEDURES 



A. ARRANGEMENT OF THE EXPERIMENT 

The measurements were made in a tank which was 300 cm 
long, 115 cm wide and 90 cm deep. The tank was filled with 
carefully de-aerated (fresh) water and had a 55 cm deep 
sand bottom (also de-aerated) . 

At first, an ITC 501 projector was used but the beam 
pattern was fairly broad leading to excessive reflections 
from the surface and the side-walls and, more importantly, 
this projector was large enough so that it would ring for a 
significant length of time compared to the desired pulse 
length. As a result, the output pulse was too long regard- 
less of the length of the pulse used as input. 

These problems led to the selection of a USRF F-41 
projector, which is a directional laboratory standard trans- 
ducer used in acoustic reciprocity calibrations over the 
frequency range from 15 to 150 kHz. This projector consists 
of lead zirconate-titanate elements cemented to a high- 
density kennametal disks. The resulting array of elements 
is about 3.8 cm wide and 5 cm high. 

The receiving hydrophone used was an omnidirectional 
LC-10 hydrophone mounted on a carriage that was free to roll 
along tracks installed on the long sides of the tank. In 
this way the horizontal distance (X) could be easily varied. 
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A block diagram of the electronics used is shown in 
Fig. 3.1. The digital oscilloscope was used to measure the 
time differences between various pulses. In order to reject 
unwanted noise, a filter was used to pass only 100 kHz to 
200 kHz. 

The choice of pulse length was governed by the criteria 
that it be (1) long enough to reach steady-state conditions 
in the projector and (2) short enough to resolve the indi- 
vidual reflections from the layers. The best compromise 
turned out to be four cycles of a 150 kHz sine wave, but, 
unfortunately, this was not short enough to completely 
define the thin plate. At 150 kHz, one cycle is 6.7 yS 
so the four-cycle pulse was about 27 yS long. A typical 
received signal using the four-cycle pulse is shown in 
Fig. 3.2 while Fig. 3.3 shows the same situation with an 
eight-cycle pulse. The four-cycle pulse gives much better 
definition of individual arrivals . 

Before the actual measurements were made, several tests 
were performed in order to insure that side-wall and surface 
reflections would not interfere with the measurements. The 
depths of the projector and hydrophone were determined by 
trial and error not only to minimize interfering reflections 
but also to optimize the separation between the desired 
arrivals. In addition, this set a practical limit on the 
horizontal range. 

An attempt was made to move the axis of the source-to- 
receiver line toward the tank wall so that two parallel 
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Figure 3.1. Block Diagram 
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Figure 3.2. 4 Cycle Pulse Signal 




Figure 3.3. 8 Cycle Pulse Signal 



measurements could be made without disturbing the sand to 
bury each plate. This introduced unacceptable side-wall 
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Figure 3.4. Location of Projector, Hydrophone 
and Plates 




Figure 3.5. Projector Is Not on the Center Line 
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reflections so, instead, three plates were buried in line 
along the long axis of the tank with the projector and 
receiver positioned as shown in Fig. 3.4. 

Here, the plates were buried at different depths and the 
projector and hydrophone were positioned as shown. (Notice 
that in the third set-up on the right the source and receiver 
are reversed to avoid reflections from the end-wall.) 

Figure 3.5 shows the side-wall reflections that are 
introduced when the source-to-receiver axis is displaced 
only 15 cm toward one side of the tank. This result should 
be comparable to Fig. 3.2 but only a position along the 
centerline of the tank turned out to be acceptable. 

B . EXPERIMENTAL PROCEDURE 

Once the experiment was set up, the first step was to 
identify the individual arrivals from a particular buried 
plate. By considering the projector depth S, the hydrophone 
depth h, the water depth H, the major incident angle, and 
the slant distance R, the time of arrival of the various 
reflections was determined by simple image theory [Ref. 3, 
pp. 427-429]. A representative diagram of the various 
arrivals is shown in Fig. 3.6 for constant sound speed in 
the water layer. 

Since a directional projector was used and was aimed 
downward toward the sediment, the surface-reflected path 
7 could be neglected. For most of the measurements, paths 1, 
3 and 4 were used. 
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Figure 3.6. Waves in the Isospeed Shallow 
Water Channel 

Basically, three methods were used to initially identify 
the arrivals on the digital oscilloscope. First, the 
arrival times predicted by the method of images were com- 
pared to the measured arrival times. The amplitude of 
the arrivals also helps in the identification process since 
the pressure amplitude is reduced by the reflection coeffi- 
cient at each of the interfaces of the layers. Finally, a 
path that reflects from the surface can be readily identi- 
fied by disturbing the water's surface with a ruler. Any 
surface-reflected arrivals change greatly on the oscillo- 
scope display and can easily be separated from those 
arrivals that only reflect from the bottom. 
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A typical display on the oscilloscope is shown in Fig. 
3.7. Here, the projector is inclined so that the beam has 
an incident angle of 55 degrees and the water depth is 35 cm. 
The arrivals shown are; 

1. direct wave, 

2. surface-reflected wave, 

3. bottom-reflected wave, 

4. wave reflected from the top of the plate, 

5. wave reflected from the bottom of the plate after 
passing through the plate, 

6. bottom-surface-reflected wave. 
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Once the arrivals are identified for one source-to-receiver 



range, the signals can be followed by eye if the receiver 
carriage is moved smoothly. 

Next, the range X is varied by moving the carriage and 

the arrival times T of the individual reflections are 

recorded. These measurements form the basis of the T-X and 
2 2 

T -X methods. When these measurements are plotted as 

arrival time T as a function of range X, the direct arrivals 

appear as straight lines and the reflected arrivals (either 

from the top of the sand layer or from any of the layers 

below) appear as parabolas. (They are actually closer to 

hyperbolas but, over the limited range available in this 

experiment, parabolas were used as an approximation.) 

2 

The T and X measurements were also plotted as T as a 
2 

function of X . On these coordinates, any of the lines are 

straight (or nearly so). From the slope of these lines, 

the sound speed either in the water or in the layers could 

be calculated. A linear regression method was used to 

compute the slope for the individual arrivals. 

Next, the ray parameter P was calculated. This value 

is fixed for a given measurement. The incident angle was 

determined from the source-receiver geometry and the T-X 

data were fit with parabolas. Then P could be calculated 

2 

directly from the parabola equation: T = aX + bX + c , 

dT/dX = P = 2aX + b. 

Once the parabolas have been determined, the layer depth 
can be found from the T-X measurements. For example, the 
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arrival from the wave that passed through the upper sand 

layer would be described by a parabola of the form 
2 

+ b^X^ + and the reflection from the top of the 

2 

sand layer would be given by T 2 = + b 2 X 2 + c 2’ The 

ray parameter P must be the same for each arrival; there- 
fore, for the wave passing through the sand layer, 

P = 2a^X^ + b^ or X^ = (P-b^)/2a^, and, for the wave reflected 
from the top of the sand layer, P = 2a 2 X 2 + b 2 or X 2 = (P-b 2 )/2a 2 
Since AX = X^-X 2 , 



P-b 2 P-b, 

AT = t {a l ( Tif> +b l ( ^I7 



P-b 9 - P-b 9 

} + c l } - {a 2 { ~2l7 ) +b 2 ( ^f )+c 2 }] 



These two values (T and X) are necessary to calculate the 

2 2 

depth and the sound speed (to compare with the T -X deter- 
mination of sound speed) of the layer. The values of layer 
depth Az and sound speed c^ are given by, 

1 22 
Az = 2 -c^AT(1 -P c^) 



2 2 

The first value of c^ is from the T -X method and AT is 
from the parabolic fit of the T-X data. The new value of 
c^ can be compared to the old value as a check on the two 
methods . 



C. SUMMARY OF EXPERIMENTAL PROCEDURE 

1. Bury the thin homogeneous plates at various depths; 

2. Identify the arrivals on the oscilloscope; 
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3 . 



Measure the horizontal distance X and the arrival 
time T of the various arrivals; 

4. Make the T-X and T^-X^ plots; 

5. Find the sound speed c, in the layer by using the 
T^-X 2 method; 

6. Determine the incident angle from the source- 
receiver geometry; 

7. Calculate the ray parameter P (fixed for each 
experiment) ; 

8. Fit the T-X data with parabolas; 

9. Calculate AX and AT by the ray parameter method; 

10. Calculate the layer depth Az; 

11. Calculate c^ and compare to the value determined 
previously. 
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IV. ANALYSIS AND RESULTS OF EXPERIMENTS 



A. MEASUREMENT AND ANALYSIS 

A measurement of the sound speed in the water was 

required for each of the reflection analyses so this was 

2 2 

measured using the T -X method. The slope of the straight 
line fit to the data points gave a sound speed of 1447 m/s 
and the correlation coefficient associated with the linear 
regression was almost 100%. For this measurement, the pro- 
jector beam was aimed horizontally and the hydrophone was 
positioned at the same depth as the projector. 

In each of the subsequent experiments, the projector and 
hydrophone depths were chosen to optimize the separation 
and clarity of the reflection signals. In every case, the 
source signal was four cycles of a 150 kHz sine wave. 

1 . Measurement 1 

For this experiment an aluminum plate (thickness of 
1.7 cm) was buried 3.5 ± 0.05 cm below the sand's surface 
and the projector depth was 8.5 cm while the hydrophone 
depth was 6.5 cm. 



WAVE 


TABLE I 

Speed of Sound (Measurement 1) 

SPEED OF SOUND (m/s) CORRELATION (%) 


ERROR (%) 


SAND 


1661 


98.3 




BOTTOM 


1415 


99.4 


2.2 


DIRECT 


1463 


99.8 


1.1 
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Table I gives the speed of sound calculated from the various 

arrivals along with the correlation coefficients of the 

linear regressions. The speed of sound in the water should 

be compared to the value of 1447 m/s discussed above. The 
2 2 

T-X and T -X data are shown in Fig. 4.1. 

As discussed in the previous chapter, the T-X data 
were fit with a parabola for simplicity. The data for the 
wave that reflected from the top of the aluminum plate was 
fit with 



T = 0.0926 X 2 - 2.89 X + 459 

and the data for the wave that reflected directly from the 
top of the sand was fit with 

T = 0.113 X 2 - 3.03 X + 410 

The geometry selected for the ray parameter determination 
gave an incident angle of 29 degrees which corresponds to a 
ray parameter P equal to 3.32 x 10 ^ s/m. This value of slope 
locates a point on each of the two parabolas. The differ- 
ence in X values gives a AX of 5.47 cm. The parabolic form 
for the reflection from the top of the plate gives the time 
difference AT of 51.4 yS. 

Next, the layer thickness Az was calculated from 
Equation (14). In this case, Az was determined to be 3.56 cm. 
This differs by only 0.3% from the measured value for the 
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overlying sand layer. (The burial depth of the plate gives 
the thickness of the top sand layer.) The value for the 
layer's sound speed was also calculated from the ray param- 
eter method but its value was not significantly different 

2 2 

from that calculated by the T -X method. 

2 . Measurement 2 

For the second measurement, the hydrophone depth was 

6.4 cm, the projector depth was 8.5 cm and the aluminum plate 

(1.7 cm thick) was buried at a depth of 2.9 ± 0.1 cm in the 

sand. The calculated sound speeds are given in Table II and 
2 2 

the T-X and T -X data are shown in Fig. 4.2. 



TABLE II 

Speed of Sound (Measurement 2) 



WAVE 


SPEED OF SOUND (m/s) 


CORRELATION (%) 


ERROR 


SAND 


1647 


94.3 




BOTTOM 


1419 


99.6 


1.96 


DIRECT 


1437 


99.1 


0.69 



Figures 4.3 and 4.4 show representative displays from 
the oscilloscope at two of the ranges at which measurements 
were taken. 

The direct wave (A) , the surface-reflected wave (B) , 
the wave that reflects from the top of the sand (C) , the 
wave that reflects from the top of the buried plate (D) , and 
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Figure 4.2. 
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(Measurement 2) 



36 



X (CM ) 





Figure 4.4. Measurement 2 at X 



39.6 cm. 



37 






the wave that reflects from the top of the sand and again 
from the water surface (E) are shown in these figures. The 
vertical line marks the beginning of the wave that is 
reflected from the top of the plate and the number gives the 
arrival time. 

In this case, the curve for the arrival reflected 
from the top of the plate is, 

T = 0.473 X 2 - 32.0 X + 1020 

and the curve for the arrival reflected from the top of the 
sand is. 



T = 0.136 X 2 - 5.72 X + 471 

_ g 

Here, the incident angle was 28.3 degrees so P was 3.30 * 10 
s/m. The ray parameter calculations gave 3.87 cm for AX and 
54 yS for dT so the layer depth Az was 3.73 cm. The actual 
sand layer thickness was 3.4 cm so the error was slightly 
less than 10%. The value for sound speed in the water by the 
ray parameter method was 1474 m/s which is about 10% high. 

3 . Measurement 3 

In this measurement, an aluminum plate 0.1 cm thick 

was buried to a depth of 2.3 ± 0.1 cm. The projector depth 

was 8.5 cm and the hydrophone depth was 5.1 cm. The result- 

2 2 

ing speed of sound measurements (determined by the T -X 
method) are summarized in Table III and the T-X data and the 
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T -X . 


data are shown in Fig. 


4.5. Figures 4.6 


and 


4.7 


show representative displays 


from the oscilloscope 


at two 


of the 


ranges at which 


measurements were taken 


o 








TABLE III 








Speed of 


Sound 


(Measurement 3) 






WAVE 


SPEED OF SOUND 


(m/s) 


CORRELATION (%) 


ERROR (%) 


SAND 


1649 




95.2 






BOTTOM 


1471 




99.9 




1.7 


DIRECT 


1510 




99.8 




4.3 



The parabola for the arrivals reflecting from the 
top of the plate was determined to be 

T = 0.290 X 2 - 13.3 X + 580 



and the parabola for the arrivals directly reflected from 
the top of the sand is 

T = 0.0484 X 2 + 0.528 X + 357 

For the optimum geoemtry , the incident angle was 27 degrees 

— 6 

so that the ray parameter P was 3.00 x 10 

The ray parameter analysis gave 2.53 cm for AX and 
33.38 ys for AT; therefore, the sand layer depth was computed 
to be 2.42 cm: less than one percent from the actual value. 
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Figure 4.5. 



T-X and T 
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Plots 



(Measurement 3) 
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Figure 4.7. Measurement 3 at X = 29.3 cm. 
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The calculated value for the sound speed in the sand layer 

was 1578 m/s and this differs by only 3 percent from the 

2 2 

value obtained by the T -X method. 

4 . Measurement 4 

For this measurement, an aluminum plate 0.3 cm thick 

was buried to a depth of 1.8 ± 0.1 cm. The projector depth 

was 8.5 cm and the hydrophone depth was 3.5 cm. The speed 

of sound measurements is summarized in Table IV and the T-X 
2 2 

data and the T -X data are shown in Fig. 4.8. Representative 
displays from the oscilloscope at two ranges are shown in 
Figures 4.9 and 4.10. 



TABLE IV 

Speed of Sound (Measurement 4) 



WAVE 


SPEED OF SOUND (m/s) 


CORRELATION (%) 


ERROR 


SAND 


1582 


99.7 




BOTTOM 


1458 


99.7 


0.7 


DIRECT 


1461 


99.9 


1.0 



The parabolic form for the reflections from the top 
of the plate is 

T = 0.0259 X 2 + 0.361 X + 429 

and the form for the reflections from the top of the sand is 
T = 0.0328 X 2 + 1.17 X + 380 
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Figure 4.8. 



T-X and T 2 -X 2 Plots 



(Measurement 4) 
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Figure 4.10. Measurement 4 at X = 47.7 cm. 
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For the optimum geometry, the incident angle was 50 degrees 

_ g 

so that the ray parameter P was 5.3 x 10 

The ray parameter analysis gave 5.80 cm for AX and 

40.4 ys for AT; therefore, the sand layer depth was computed 

to be 1.74 cm. The calculated value for the sound speed in 

the sand layer was 1647 m/s. In both cases, the results 

compare very well with the measured value for layer depth 
2 2 

and the T -X value for sound speed. 

5 . Measurement 5 

For this measurement, the projector depth was 8.5 cm 

and the hydrophone depth was 3.5 cm. The aluminum plate 

(0.1 cm thick) was buried at a depth of 5.0 ± 0.1 cm. The 

speed of sound measurements are summarized in Table V and 

2 2 

the T-X data and the T -X data are shown in Fig. 4.11. 
Samples of the received arrivals at two ranges are shown in 
Figures 4.12 and 4.13. 



TABLE V 

Speed of Sound (Measurement 5) 



WAVE 


SPEED OF SOUND (m/s) 


CORRELATION (%) 


ERROR 


SAND 


1577 


99.9 




BOTTOM 


1468 


100 


• 

1— 1 


DIRECT 


1488 


100 


00 

• 

CM 



The parabolic form for the reflections from the top 
of the plate is 
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Figure 4.11. 



T-X and 



X 



Plots 



(Measurement 5) 
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Figure 4.12. Measurement 5 at X = 28.2 cm. 




Figure 4.13. Measurement 5 at X = 48.8 cm. 
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T = 0.0450 X 2 - 0.213 X + 407 

and the form for the reflections from the top of the sand 
is 



T = 0.0463 X 2 + 0.562 X + 339 

The incident angle for determining the ray parameter was 

— fi 

33.4 degrees so that the ray parameter P was 3.70 x 10 

The ray parameter analysis gave 9.59 cm for AX and 

71.5 ys for AT; therefore, the sand layer depth was com- 
puted to be 4.58 cm — about six percent below the actual 
value. The calculated value for the sound speed in the sand 
layer was 1904 m/s which was significantly higher than any 
other of the determinations . 

B. OVERALL RESULTS 

In all, seven experiments were performed but the data 
from two of the experiments was not good enough to analyze. 
The other five experiments have been reported above and, 
except for number 5, the results were quite consistent with 
reasonable errors. 

Table VI gives a summary of the results of the five 

experiments that were analyzed and compares the results of 
2 2 

the T -X method and the ray parameter method. The average 
error in determining the depth of the layer was 3.5 percent 
so this measurement can be done effectively; however, the 



48 



